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Abstract 

We revisit the QCD sum-rule treatment of the isospin-breaking correlator 
(0|T[^(x)^(0)]|0), in light of the recent claim that a previous treatment 
produced results incompatible with known chiral constraints. The source of 
the error in the previous analysis is identified, and a corrected version of the 
sum-rule treatment obtained. It is then shown that, using input from chiral 
perturbation theory, one may use the resulting sum rule to extract information 
on the leading chiral behavior of isospin-breaking parameters associated with 
the coupling of excited pseudoscalar resonances to the axial currents. A rather 
accurate extraction is possible for the case of the r]'. Demanding stability of 
the sum-rule analysis also allows us to improve the upper bound on the fourth- 
order low-energy constant, Lj. 
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I. INTRODUCTION 



One of the most attractive features of chiral perturbation theory (ChPT) [|l|] is that it 
provides a framework for constructing effective hadronic Lagrangians in the most general 
possible way that implements both the symmetries of QCD and the symmetry breaking 
pattern of the approximate chiral symmetries of QCD. As such, it fully incorporates the 
consequences of QCD in the low-energy regime. The price to be paid for using only symmetry 
arguments is that every term in the effective Lagrangian, C e s, allowed by these arguments 
will appear, multiplied by an undetermined constant (referred to as a low-energy constant, 
or LEC). These LEC's could, in principle, be computed from QCD, but must be treated as 
parameters to be determined phenomenologically if one does not go beyond ChPT. 

Although such effective Lagrangians are necessarily non-renormalizable, Weinberg's 
counting argument 0j shows that only a finite number of terms in C c r contribute if one 
expands to fixed "chiral" order, that is in powers of the external momenta (generically de- 
noted as p) and current quark masses (where m q counts as order 0(p 2 )). As a result, in the 
chiral expansion of any low-energy observable, the general form of the dependence on exter- 
nal momenta and light quark masses, to a given order, can be computed straightforwardly 
from the form of the relevant terms in £ e g. 

Since this formal dependence is a rigorous consequence of the symmetries and approxi- 
mate symmetries of QCD, it follows that ChPT can be used to place constraints on treat- 
ments of the same observable using other methods. Indeed, if one makes a chiral expansion 
of the results obtained by any other method and finds that terms present in ChPT to a 
given order are missing, then one knows unambiguously that either the method itself, or 
some truncation employed in it, is incompatible with QCD. This is true regardless of the 
rapidity of convergence of the chiral series in question: all terms required by the symmetries 
of QCD must be present if the method is to correctly incorporate the consequences of QCD. 

An example of the use of such constraints is provided by the analysis of the nucleon mass 
using QCD sum rules. Standard treatments were shown to produce an expression for in 
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terms of condensates that implies the presence of certain chiral logarithms inmjv, although 
such contributions are known from ChPT to be absent The source of this problem was 
found to be a failure to treat properly the contribution of the ttN continuum to the sum rule 
in the original analyses [|J ; including the leading contributions from such states restores the 
correct chiral behavior of m^. 

A more severe problem of the same type has been pointed out in the case of the isospin- 
breaking axial correlator 

nj(ff) = ijd A x e^(0|T(4(x)^(0))|0> = n^ 2 )^ - n 2 ( g 2 )gV , (i) 

where A 3 , A s u are the neutral members of the octet of axial currents A ^ = 97^75 ^q. This 

correlator was first analyzed using QCD sum rules in Ref. (CHM). As shown in Ref. ||, 

however, if one writes Il 1 (g 2 ), which contains the 7r° and rj pole contributions, in the form 

9v 9* \ | = / q 2 (9y - 9n) + {g^m 2 - ^m 2 ) \ | 

q 2 — m 2 q 2 — m 2 J \ (q 2 — m 2 )(q 2 — m 2 ) J 

then the expression for g v — g n (given by the slope of the numerator with respect to q 2 ) 
obtained from the sum-rule analysis is lacking both the leading analytic and leading non- 
analytic terms from its chiral expansion in terms of the light quark masses. (The demon- 
stration of this is reviewed briefly below in Sec. II.) 

In this paper we revisit the sum-rule analysis of the axial correlator above, and identify 
the source of this problem. We then obtain a corrected version of the relevant sum rule and 
show how it can be used to extract information on isospin-breaking couplings of the higher 
pseudoscalar resonances. 

The paper is organized as follows. In Sec. II, we briefly review the sum-rule and ChPT 
analyses of the correlator. In Sec. Ill, we identify the problem with the previous sum-rule 
treatment and work out the corrected version of the relevant sum rules. In Sec. IV, we 
show how one can use information from ChPT as input into the sum rule. We also clarify 
the physical content of the corrected sum rule, extracting in the process information on the 
isospin-breaking couplings of the higher pseudoscalar resonances to the axial currents. We 
conclude in Sec. V with a brief summary. 
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II. PREVIOUS CHPT AND SUM-RULE TREATMENTS 



We provide here only a very brief review, which will serve also to fix notation. For more 
details the reader is referred to Refs. [f| and || for the sum-rule and ChPT treatments, 
respectively. 

We first review the sum-rule treatment. [] As usual, the aim is to write dispersion relations 
for n^g 2 ) and n 2 (g 2 ) which relate integrals over the relevant physical spectral functions to 
the behaviors at large spacelike q 2 , where the operator product expansion (OPE) becomes 
valid. One then Borel transforms the resulting dispersion relations in order to exponentially 
suppress the higher-energy portions of the spectral integral on the phenomenological side 
and simultaneously factorially suppress the contributions of higher dimension operators on 
the OPE side. The scalar correlators Ui(q 2 ) and n 2 (g 2 ) in Eq. (|1|) have been chosen in such 
a way that, from the asymptotic behavior of Il^(g) in QCD, it is known that the relevant 
spectral integrals converge without subtraction. Note that the definition of II2 employed 
here agrees with that used in Ref. ||, but differs from that in Ref. || by a factor of — q 2 . 

On the phenomenological side, the axial- vector resonances contribute to both Ili(g 2 ) and 
n^g 2 ). In the narrow- width approximation, their contributions to the complete spectral 
function are written 



The pseudoscalar resonances, in contrast, contribute only to IL;. Following the convention 
of earlier works, we write these contributions as 



-(imn^g 2 ))^ = gj(q 2 - m 2 ) - g v 5(q 2 - m 2 ) + g n >5(q 2 - m 2 ,) + g^5{q 2 - m 2 ,) + ■ ■ • . 



(The minus sign in front of g v is conventional and related to the fact that, so defined, g v = g 
at leading order in the chiral expansion.) 




(3) 



(4) 



See, for exampL 



e, Refs. 0-10 1 for details of the general method of QCD sum rules. 
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On the OPE side of the sum rule, the expressions for the scalar correlators have been 
worked out by CHM up to and including operators of dimension 6, and to order 0(m q , a s ). 
Neglecting electromagnetic effects, the results have the form (with Q 2 = —q 2 ) 



ni(g 2 ) 
n 2 (g 2 ) 



1 



4x/3 
1 

4^3 



C lnQ 



c. 



^ q 2 ^ g 4 ^ g 6 

Going + _-_ + _ 



(5) 



where Cq and C 2 vanish at the level of the truncations noted above, and 



C 2 = 2[m u (uu) - m d (dd)] , 
C 3 = ^tto> s [(uu) 2 - (dd) 2 } 



(6) 



If one were to include higher-order terms in Eq. (|^), Co would receive contributions at order 
OiptEM-, 0i si iriq) and C\ at 0(m 2 ). An argument analogous to that of Shifman, Vainshtein 
and Zakharov || for the corresponding isospin-conserving correlator H^ u (q) shows that the 
higher-dimension operators not included in these expressions are also all explicitly of order 
0(m 2 ). The form of the dimension 6 coefficient C3 in Eq. @ has been obtained assuming 
vacuum saturation. 

As can be seen from the Lorentz structure of Eq. it is possible to remove the contri- 
butions of the axial- vector mesons by considering the combination 

np^n^-n^ 2 ) . (7) 

CHM, motivated by this observation, write a dispersion relation for Up(q 2 ) in the form 

1 ImIIp(s; 



U P (q 2 



ds . 



7T s — q* 

When Borel transformed, this relation gives CHM's sum rule, 



(8) 



2C, 



9lx e- m "l M2 - g v e""^ 2 



+ 



(9) 



ViWV M 2 

where M is the Borel mass parameter and the dots refer to the contributions of higher 
pseudoscalar resonances. CHM then neglect higher resonance contributions and use this 



sum rule, together with its derivative with respect to the Borel mass, M, to solve for g n and 
g n . This procedure leads to their result 



{g v - sw) 



CHM 



m 2 — ml 



(10) 



,2^3M 2 y 

for the slope of the numerator on the RHS of Eq. (Q). 

The analysis of U 3 ^ u (q) at next-to-leading (1-loop) order in ChPT is straightforward, and 
follows standard methods. We employ throughout the notation of Gasser and Leutwyler [|l| . 
The result for IT 2 (g 2 ), recast so as to correspond to the definition employed in this paper, is 



j-j-l— loop ^2 



B (m d - m u ) 



32tt 5 



(log(mLV) + l)-8L^ 2 



(11) 



where B is the usual second-order LEC, related to the quark condensate in the chiral limit, 
fi is the renormalization scale, and L r 5 (fi 2 ) is a renormalized fourth-order LEC. Note that 
n 2 (g 2 ) results solely from contact terms (that is, terms in £ e g that are quadratic in the 
external axial sources). To this order, Ili(g 2 ) is saturated by the tt° and 7] pole terms. From 
a similar analysis, one finds, for the coefficients g n and g n appearing in Eq. (|J), 



9n = flt\ and g v = / 2 e 2 



(12) 



where F n , f v are the physical n, i] decay constants and ei, e 2 are isospin-breaking parameters 
defined by 



(0|4|tt) = iUe iqit 
(0\Al\ V ) = -if v e 2qfl . 



(13) 



The expressions for /„., f v , e\ and e 2 valid to 1-loop order can be found in Ref. |1[. 

The problem with the sum-rule treatment is exposed when one uses the known chiral ex- 
pansions of the meson masses and quark condensates to rewrite the sum-rule result, Eq. ([T0"D. 

as 



{% ~ 9*) 



CHM 



9 <> F 9 



8 Bq(7ti s — m)(m s + 2m) 



M 4 



+ 



(14) 



to order 0(m 2 ). Here F is a second-order LEC, equal to f w in the chiral limit, and 8q is the 
leading-order 7r°-7/ mixing angle, 

4 V m s — m J 

with rh = (m u + m^)/2. Comparing this expression with the corresponding one obtained 
from the 1-loop ChPT results, 

/ x r, ^9 f ( m l — ™hc) , / 9 ; 9x B (m s — rh) 

3250(771, -m) , 2 n \ ( , a \ 
+ ^ L 5 (fi ) + ■■■], (16) 

one sees that the sum-rule expression is lacking both the leading analytic and leading non- 
analytic terms in its chiral expansion 0, and hence is incorrect. Moreover, the numerical 
consequences of this are significant: the sum-rule value for the slope is more than an order 
of magnitude smaller than that given by ChPT. 



III. CORRECTED VERSION OF THE SUM-RULE ANALYSIS 

The key to understanding the origin of the problem with the CHM sum-rule analysis lies 
in Eq. (§). This relation follows from general properties of analyticity and unitarity under 
two assumptions: (a) that the singularities of ilp(g 2 ) consist solely of those associated with 
physical intermediate states and (b) that Ilp(g 2 ) converges sufficiently fast that no subtrac- 
tions are required. The latter assumption is explicitly verified by the known asymptotic 
behavior of n^g 2 ) and n 2 (g 2 ) in QCD. The former, however, is more subtle, since there 
can also be singularities of purely kinematic origin. In the case at hand, Eq. (|TT| ) shows 
explicitly that n.2(q 2 ) has a kinematic pole at q 2 = 0. As a consequence, the correct version 
of the dispersion relation Eq. (|8p must include the contribution of this kinematic pole to the 
underlying contour integral. Another way of saying this is that it is q 2 U2{q 2 ) which satisfies 
a dispersion relation without kinematic pole terms. The dispersion relation for this function, 
however, requires one subtraction in order to converge. The resulting subtraction constant 
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gives rise to the kinematic pole term of I^g 2 ). Its value is calculable in ChPT, and turns 
out to correspond precisely to the contact contributions given in Eq. (|TT|). 

Bearing this in mind, it is straightforward to write down the corrected dispersion relation 
for U P (q 2 ), 



2 1 B (m d - m u ) 
LLp{q ) - 



q 2 



32 



. 1 Imllp(s) , . . 

+ / ^rds , 17 

' 7T s — q 2 



where Imilp(s) includes only the spectral strength associated with pseudoscalar states. The 
corresponding Borel-transformed sum rule is then 



Co 



2V3M 2 



B (m d - m u )\ ( 3 



,2 /» /f 2 



/M* 



- 9v e-< IM '+9, e-^ + g* e~<>^ + 



(18) 



As one might expect, the inclusion of the kinematic-pole contribution cures the problem 
of the incorrect chiral behavior of g v — g-n- To see this, consider the O(M ) terms of Eq. (|T8|). 
Bearing in mind that g v i, g w i, ... are all of order 0(m 2 ) ||, one has 

B (m d - m u )\ ( 3 



= [~ ^ U, j (3^ (log«//i 2 ) + l) - 8Ll(^)) +g„-g v + 0{m\) , (19) 



where the first term on the RHS results from the kinematic pole in Eq. (1171) . Without this 
term, one gets g v — g n = 0(m 2 ), as found by CHM. In contrast, using the corrected sum 
rule, one finds that Eq. (|T9|) is simply an alternate form of Eq. (|16|) , as required. 

To clarify the physical content of the remaining pieces of the sum rule, Eq. ([T8|), it is 
useful to note the chiral order of various quantities appearing therein. In particular, the chiral 
expansions of g n , g^, m 2 , and m 2 , start at order O(p ), m 2 , m 2 and C 2 at 0(p 2 ), and (as 
already noted above) g v >, g n i at 0{m 2 q ) = 0{p 4 ). After the cancellation embodied in Eq. (|19D, 
the only 0(p 2 ) terms remaining in Eq. flTB| ) are those in C2 and —g^m 2 + g v fn 2 - Using the 
leading-order expressions (uu) = (dd) = —B F 2 and g n = g v = 9 F 2 , it is straightforward to 
show that the 0(p 2 ) terms on both sides of the sum rule also match properly. To this order, 
the matching is just an isospin-breaking version of the Gell-Mann-Oakes-Renner relation. 



S 



The information obtained in the previous paragraph is all that we can extract from 
Eq. (|I5D in its present form. This is because terms of 0(p A ) have not been included on the 
OPE side of the sum rule. If one wishes to use the sum rule to obtain information about 
anything beyond the leading and next-to- leading order behavior of g n and g v , one must, 
therefore, restore the 0(m 2 ) terms to the OPE. This is easily accomplished starting from 
the expression for the corresponding terms in the OPE of the analogous isospin-conserving 
correlator, as given in Ref. JlTJ. The result is 



n P (g 2 



OPE 



V3Q 2 



3(m 2 d 



7T 2 V yU 2 



~^ rrid(dd) — m u (uu) 



Q 2 



{m 2 d - m 2 u ) 2 
H — \ot s G ) 



TlQ 4 



(20) 



where we have kept terms only up to dimension 4 and written down the coefficient func- 
tions only to leading order in a s . Substituting the expression [np(g 2 )] OPE into the LHS of 
Eq. ([T7|) and Borel transforming, we obtain an improved version of the corrected CHM sum 
rule, Eq. fll^p. To facilitate subsequent analysis, it is convenient to multiply both expressions 
for the correlator by Q 2 before Borel transforming (thereby eliminating the contribution of 
the kinematic pole). We also follow standard practice and introduce a continuum threshold 
parameter, sq, representing the point beyond which the hadronic spectral function is mod- 
elled by its perturbative QCD counterpart. The contribution corresponding to the integral 
over that portion of the phenomenological spectral function can then be moved to the OPE 
side of the sum rule. The result of these manipulations is the sum rule, 



2 , ' /J - n n, 2 , -</M* + gpm 2 p e -"4/M 2 



i 2 IM 2 

g^m-^ e - g v m- e 



1 

7! 



^|^|) M2 ( e -WM 2 _ x) + f m _ mu{uu)) _ ^fl {as G 2 ) 



(21) 



where the sum on the LHS now runs over pseudoscalar resonances with squared masses less 
than s - 

The chiral expansion of the sum rule, Eq. (pi]), contains terms of order 0(p 2 ) and higher, 
together with the usual chiral logs, which start at order 0(p 4 \np). The 0(p 2 ) terms are the 
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same as those in Eq. (|1|) and so it is easy to see the sum rule is consistent to this order. 
Since only the light-quark condensate and the quantities g n , g v , m 2 and m 2 contain leading 
chiral logs, these contributions must also cancel in Eq. (|2"ID (as verified below). Finally, the 
expansion of g^rn^ e~ m ^l M2 — g v m 2 e~ m ^ M2 to order 0(p 4 ) can be found from known 1-loop 
expansions, and that for m^dd) — m u (uu) can be obtained from a straightforward 1-loop 
calculation. With these results, we may employ this sum rule to obtain a relation describing 
the leading chiral behavior (0(m 2 ) = 0(p A )) of the isospin-breaking parameters g v >, g n >, . . . 



for the heavy pseudoscalar mesons. No further information can be extracted from Eq. (|2lD 
without 2-loop ChPT calculations as input. 

To verify the cancellation of the chiral logs, and to obtain the promised sum rule for the 
leading chiral behavior of gy, gy, . . . , we expand the tt, rj and condensate terms in Eq. fl21~|) 
to order 0(p A ). To do so for the 7r and rj terms appearing on the LHS requires only the 
1-loop expressions for f n , f v , e\ and €2 given by Gasser and Leutwyler |L|. The results are 



g n m 2 n - g v m 2 v 



-e F 2 



B (m s 



128B 2 (m 



m ) - B ( 

rh) r 



B\m 2 
6n 2 F 2 



3F 2 

log (m 2 K /fi 



'9{m s — 2rh)£ n + 6m s £ K + {m s + 2m)(i ri \ 
x 72tt 2 F 2 J 

m s + 2m)L r 6 (fi 2 ) - 2{m s - m)U 7 + 2mL r s (fi 2 ) 
2 \ B 2 m(m s — fn) 



12n 2 F 2 



(22) 



and 



-p (-9nmi + 9 v ml 



M 2 



1Q6 F 2 2 2 

—BAm* — m) 

9M 2 oV s ' 



(23) 



where lp = m 2 p log (mp/n 2 ) and all other notation is as in Gasser and Leutwyler [Q. For 
the condensate contributions on the RHS, we require the expressions for (dd) and (uu) valid 
to order 0(rrid, — m u ). These are easily obtained, and can be written 



(uu) = (uu)i + 5 
(dd) = (uu)i — 5 



(24) 



where (uu)j is the 1-loop expression for the condensate in the isospin-symmetric limit, also 
to be found in Ref. [fill, and 
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lm d - m u 



Q4ir 2 F 2 (m s — rh) 



+ 



B (l + \og(m 2 K /n 2 )) 45, 



32n 2 F 2 



(25) 



From Eqs. (|2~3| ) and 

rrid(dd) — m u (uu) 



) it follows that, to order 0(p 4 

-e F 2 



m) 



B 2 rh(m s — rh) t B 2 m 2 ^ ( m 



-Bn 



12n 2 F 2 Qn 2 F 2 
9(m s — 2rh)£ n + 6m s iK + (jn s + 2rh)E ri \ 



72n 2 F 2 



+ 



GABnirrie — rh) r 



3F 2 



2(m s + 2m)Ll{tf) + 2mL r 8 (^) + mH^/j 2 " 



• (26) 



To obtain a sum rule for the leading chiral behavior of the higher pseudoscalar resonances, 
we make use of Eqs. ([22] ) and (|23|) to replace the leading terms of the 7r and rj contributions 
in Eq. (pTj) . The terms of higher order in m n and may be neglected since they are at 
least of order 0(p 6 ) in the chiral expansion, and they are numerically small for the Borel 
masses of interest. Finally, inserting the chiral expansion of the quark condensates, Eq. (pB"D, 
into this sum rule, we get 

x/3M 2 



£ 9pm 2 P e~ m r /M2 

P^TV,T) 



(e" so 



jM 



8tt 2 

64S 2 (m s - m)9 



mj + 



8V3tcM 2 

A(m s - rh)L r 7 - 2rhL r 8 (jj 2 ) + rhH^fj 2 ] 



16 



9M 2 



9 F 2 B 2 (m s 



m) 



(27) 



Note that all of the chiral logarithms have cancelled, leaving only terms that start at order 

£V). 

It is worth noting that the term involving the chiral LEC's makes a numerically signif- 
icant contribution to the sum rule and is dominated by L 7 . Moreover, phenomenological 
treatments that use resonance exchanges to generate the LEC's in the effective Lagrangian 
of ChPT PD|J21| show that L 7 receives contributions only from flavor- singlet pseudoscalar 



states. Hence it already follows, without any more detailed analysis, that the sum rule im- 
plies the existence of significant isospin-breaking rf coupling. Before presenting the results 
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of our analysis for the isospin-breaking parameters, gp, it is also worth stressing a number 
of features of the sum rule, Eq. (p7[) , which imply that, once U 7 is fixed, these results for at 
least the rj coupling should be quite reliable. 

The renormalized LEC's L 7 , L r s have been determined phenomenologically and are rea- 
sonably well-known (see, for example, Refs. [0,0 for recently updated values). The re- 
maining LEC, can be related, for example, to the isospin-breaking condensate ratio 
7 = [(dd) — (uu)}/ (uu) . This ratio has been estimated in a number of sum-rule analyses 



jl^ H^TOIlTtP^I ■ Using any of the values of 7 obtained in these treatments to estimate H^-, 
the resulting values are such that the LEC combination in Eq. (^) is dominated by the L r 7 
term. In particular, the uncertainty in the LEC combination associated with the sum of the 
errors on the phenomenological determinations of Lg and H"^ (where the latter error is taken 
to correspond to the entire range of values cited above) is an order of magnitude smaller 
than that associated with the error on the existing phenomenological determination of L r 7 . 
We may, therefore, ignore the effect of the uncertainties in the values of Lg and H%. This 
feature of the analysis results from the fact that the coefficients of Lg and Lf£ are suppressed 



by a factor of (m s — m)/m ~ 23 |19| relative to that of L r 7 . The uncertainty in the ratio 
m s /rh which enters this suppression is, of course, also completely negligible. Note that we 
do not require an explicit input value for m s since, to the order considered in the chiral 
expansion, we may take 

B 2 (m s - rh) 2 = {m 2 K - m 2 n ) 2 . (28) 

On the phenomenological side of the sum rule, we expect contributions from all of the 
higher pseudoscalar resonances, 7/(958), 77(1295), 7r'(1300), 7/(1440), 7r'(1800), .... The 
7r'(1300) is relatively broad (r = 325 MeV |22j) and spans the region between the 7/(1295) 
and the 7/(1440). Therefore, without keeping terms of yet higher dimension in the OPE, 
we have too little information in the sum rule to both adequately parametrize the spectral 
function in the region between ~ 1300 and ~ 1450 MeV and at the same time to use the sum 
rule to extract the values of all such parameters. Hence we concentrate on the extraction 
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of g V ', parametrizing the 7/(1295), 7r'(1300), 77(1440) region in terms of a single effective 
contribution of zero width located at around 1375 MeV. By varying the position of this 
contribution between 1300 and 1450 MeV, we have verified that the extracted value of g v i 
is not sensitive to this approximation, varying by ~ ±6% over this range. This is a factor 
of 6 smaller than the variation induced by the uncertainty in the input value of L r 7l which 
we discuss in more detail below. 

The effective strength parameter describing the 7/(1295), 7r'(1300) and 7/(1440) region 
(which we denote by gy in what follows) is, of course, much more sensitive to the assumed 
position of this strength. The corresponding uncertainty in the extraction of g n > is ~ 15%, 
which is significant, although still much less than the ~ 60% associated with L T 7 . The 
stability of the determination of is attributable, to a large extent, to the fact that the 
residual term proportional to L r 7 provides the major contribution to the sum rule; as already 
noted above, L r 7 is known to receive contributions only from flavor- singlet states, of which 
the 7/' is nearest and hence should provide the dominant contribution. This feature of the 
sum rule is also responsible for the greater sensitivity of g n i to the input value chosen for 
the location of the effective strength describing the 7/(1295), 7r'(1300), and 7/(1440) region: 
the combined effective contribution to the sum rule is small relative to the dominant rj term 
and the extracted value can therefore depend sensitively on the assumed separation from 
the t/' peak. 

Having employed information from ChPT to fix the low-lying 7r and 7/ contributions to 
the original sum rule, and explicitly modelled the contributions up to 1.44 GeV, we note that 
there is now a significant gap to the next resonance contribution at 1.8 GeV. We therefore 
expect that Borel masses of order 1 — 1.5 GeV will suppress the contributions of higher 
resonance on the phenomenological side of the sum rule. 

On the OPE side it turns out that the situation is also rather favorable. First, the 
gluon condensate term turns out to be numerically very small compared to the dominant 
Ly contribution. Indeed, if we take for definiteness the value for this condensate advocated 
in Ref. |23| (which is similar to that employed, for example, in Ref. [^4fl), 
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ip±G 2 ) = 0.03 ± 0.015 GeV 4 , (29) 

7T 

which includes rather conservative errors, then we find that this uncertainty corresponds to 
< 0.3% variations in g v > and g^. 

The perturbative contribution (the first term on the RHS of Eq. (p7|)) is similarly small. 
This is fortunate since recent analyses [[25-27] suggest that conventional sum- rule determi- 



nations of the light current quark masses [23,28], 29] may have overestimated these masses 
by as much as a factor of 2. For the central value of L 7 , allowing m u + to vary between 



the conventional value, 12 MeV [24] and 6 MeV produces a variation of only 2.5% in g v i 
and g n i. Such an uncertainty is again much smaller than that arising from the errors on L 7 , 
and hence can be neglected. The smallness of this perturbative contribution also implies 
that the analysis should be rather insensitive to the continuum threshold parameter, sq. We 
expect that this should lie somewhere in the vicinity of the onset of the 7r'(1800) resonance. 
In our analysis, we find, for example, that varying sq by ±1 GeV 2 about a central value 
s = 3 GeV 2 produces variations of < 1% in g v > and g^>. 

From the above discussion, we see that the RHS of the sum rule in Eq. (^7j) is dominated 
by the terms that are directly calculable using ChPT. The first of these, involving the 0(p 4 ) 
LEC's, is the piece of the quark condensate term from the OPE that remains after cancel- 
lation against ir and 77 contributions from the phenomenological side of the sum rule. The 
second consists of the remaining (9(m 2 ) tt and rj contributions from the phenomenological 
side. Numerically it is more than a factor of 2 smaller than the LEC term, for M > 1 GeV 2 , 
and of the same sign. The major uncertainty in the values of these terms is that arising 



from the phenomenological determination of the (scale-independent) LEC fl2| , |13 



L r 7 = (—0.4 ± 0.15) x 10 . (30) 

For completeness we list below the remaining input values (apart from well-determined 
meson masses): 

m u + rrid = 9 MeV 
14 



IG 2 ) = 0.03 GeV 4 
s = 3.0 GeV 
L r 8 (m 2 p ) = 0.9 x 10~ 3 
H r 2 {m 2 p ) = -7.5 x 10" 4 

77V = 1-375 ± 0.075 GeV 
m s /rh = 24.4 , 

r = md - mu = 0.3 ± 0.05 , (31) 
m d + m u 

where by ttv we mean the location of the effective strength for the 7/(1295), 7r'(1300), 7/(1440) 
region, as discussed above. In most cases we have not shown the corresponding uncertainties, 
since, as already noted, the variations in the results associated with them are small. Apart 
from Lj, the largest uncertainty is that associated with the choice ttv, which parametrizes 
the strength lying above the t/. 

Also significant is the uncertainty associated with the isospin-breaking mass ratio, r 
fl9|| . The quoted range covers a wide range of possibilities for the degree of breaking of 



Dashen's theorem [30] for the electromagnetic contribution to the kaon mass splitting. The 



recent results of Refs. |3T|-p3| would appear to confirm a larger value for the breaking, as 
suggested by earlier analyses P4}-^1, and hence larger values of r in the quoted range, with 
a somewhat smaller resulting error. Since the subject is not yet fully resolved (see Ref. ET 



for a detailed list of recent work on the subject, including some work advocating smaller 
violations of Dashen's theorem [37H), we have refrained from attempting to make a revised 
estimate for the input central value and error on r. In any case, every term on the RHS of 
Eq. (|27|) contains one factor of rrid — m u , so that this uncertainty enters only into the overall 
normalization of the final results. It does not, therefore, affect the stability analysis of the 
sum rule, and it can be removed by quoting results in the form gp/O^F 2 . 

For a given set of values for the input parameters L r 7 and 77V > we look for values of 
and g n ' that bring the two sides of the sum rule into agreement over a range of Borel mass 
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values. A convenient way to do this is to use the sum rule, Eq. (|27|), and its derivative 
with respect to M, at a fixed value of the Borel mass, as simultaneous linear equations for 
g v r and g n i. If a region is found where the results of this procedure are independent of M, 
then this indicates the existence of a stability window where the two sides of the sum rule 
match. In Fig. 1 we show some typical results for a case where we obtain good stability, 
L 7 = -0.34 x and = 1375 MeV, with g v > = 2.88 x 1(T 5 and = -5.57 x 1(T 6 . The 
two curves are essentially indistinguishable, except at the very lower end of Borel masses 
displayed. 

As \L r 7 \ is decreased, the stability window moves to larger values of M. In this region, 
the perturbatively modelled continuum becomes increasingly important in the spectral rep- 
resentation of the correlator and so the sum rule becomes unreliable for the determination of 
resonance properties. In contrast, as \L 7 \ is increased the stability window moves to smaller 
values of M and also becomes very much narrower. In fact, for values of \L r 7 \ that are larger 
than about 0.48 x 10 -3 we are unable to find a stable matching between the two sides of 
the sum rule. This occurs before the window reaches sufficiently small values of M that the 
convergence of the OPE becomes questionable. We are thus able to use the sum rule to 
make a somewhat improved determination of the LEC L 7 , reducing by about a factor of 2 
the distance to the upper bound on its magnitude compared to the ChPT result, Eq. ( |30"1) 



12,13 



For values of L 7 in the range —0.25 x 10 3 to —0.48 x 10 3 , we obtain 

g v ,/9 F 2 = 0.42 ±0.15 

g^/OoF 2 = -0.13 ±0.07 . (32) 

The dependence on r has been scaled out of these results, as discussed above, and so the 
dominant uncertainties quoted in Eqs. (B^) are those associated with the range of values for 



L 7 . Allowing for the uncertainty in r taken from [I9|, our values for the isospin-breaking 
parameters are 

g v > = (3.6 ± 1.9) x 10~ 5 GeV 2 
16 



g v , = (-1.1 ±0.8) x 1(T 5 GeV 2 



(33) 



IV. SUMMARY 

In this paper, we have revisited the sum-rule treatment for the isospin-breaking axial 
correlator, correcting the error in a previous treatment which led to the incorrect chiral 
behavior of the slope parameter g v — g w . Including the kinematic pole omitted from the 
previous treatment restores the correct chiral behavior of the correlator. We have then 
used the explicit evaluation of the n and r\ contributions to the correlator at next-to-leading 
order in ChPT to obtain a rather well-behaved sum rule for the leading chiral behavior of 
the isospin-breaking parameters, gp, of the higher pseudoscalar resonances. This sum rule 
has been analyzed and shown to provide a rather reliable estimate for g v r, once one has fixed 
the chiral LEC, L r 7 . The requirement of the stability of this sum rule is shown, moreover, to 
provide a somewhat improved determination this LEC by reducing the upper bound on its 
magnitude. 
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FIG. 1. The OPE versus the phenomenological side of the improved sum rule as a function of 
the Borel mass, M. The dotted line is the OPE side, the dash-dotted line the phenomenological 
side in units of GeV 2 . 
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